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ABSTRACT 

Magnetic field amplification in astrophysics ultimately requires an 
understanding of magnetohydrodynamic turbulence. Kinetic helicity has 
long been known to be important for large scale field growth in forced 
MHD turbulence, and has been recently demonstrated numerically to be 
asymptotically consistent with slow mean field dynamo action in a periodic box. 
Here we show numerically that the magnetic spectrum at and below the forcing 
scale is also strongly influenced by kinetic helicity. We identify a critical value, 
fh,crit above which the magnetic spectrum develops maxima at wavenumber = 1 
scale and at the forcing scale, For / < fh,crit the field peaks only at the resistive 
scale. Kinetic helicity may thus be important not only for generating a large 
scale field, but also for establishing observed peaks in magnetic spectra at the 
forcing scale. The turbulent Galactic disk provides an example where both large 
scale (> supernova forcing scale) fields and small scale (< forcing scale, with 
peak at forcing scale) fields are observed. We discuss this, and the potential 
application to the protogalaxy, but also emphasize the limitations in applying 
our results to these systems. 

Subject Headings: MHD; turbulence; ISM: magnetic fields; galaxies: 
magnetic fields; stars: magnetic fields; methods: numerical 
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1. Introduction 

The origin of magnetic fields and the dynamics of 3-D magnet ohydro dynamic (MHD) 
turbulence in astrophysical sources are problems of long standing interest (e.g. Cowling 
1957; Moffatt 1978; Parker 1979; Krause & Radler 1980; Zeldovich et al. 1984). The 
standard in situ mean field dynamo (MFD) model (Moffatt 1978; Parker 1979; Krause & 
Radler 1980) of the large scale (= scales greater than the turbulent forcing) magnetic field 
origin can be thought of as a framework for understanding an inverse cascade of magnetic 
helicity, initiated by a forcing of kinetic helicity (Pouquet et al. 1975). While the role of 
kinetic and magnetic helicities are important for in situ non-local inverse cascade models 
of large scale fields, or MFDs, the small scale (= scales at or below the turbulent forcing 
scale) dynamo does not explicitly require helicity to amplify total magnetic energy density 
(Zeldovich et al 1983; Parker 1979). Non-helically forced turbulent amplification of the 
small scale fields, and fully helical forced growth of large scale fields have been recently 
simulated (Cho & Vishniac 2001; Chou 2001; Brandenburg 2001; Maron & Cowley 2001). 

But there is an important subtlety which has not yet been addressed. Though 
numerical work generically shows that the total energy of the small scale field in turbulent 
media saturates to of order equipartition with the kinetic energy spectrum, non-helical 
small scale dynamos produce a peak of the magnetic energy spectrum on the resistive scale 
for magnetic Prandtl number > 1, not on the forcing scale (Chou 2001; Maron & Cowley 
2001). This contradicts, for example, observations of the Galactic magnetic field which has 
a peak in the spectrum on the turbulent forcing scale (Beck et al. & 1996) and maximally 
helical simulations (Brandenburg 2001). Here we show that forcing with varying levels of 
fractional kinetic helicity affects the overall spectral shape at large and small scales. 

In section 2 we discuss the equations and the simulations. In section 3 we give the 
results and in section 4 the interpretations. We conclude in section 5. 

2. Equations and Numerical Scheme 

We investigate forced helical MHD turbulence. We write the magnetic field in velocity 
units and so define b = B/V^r, where B is the magnetic field. Incompressibility is assumed 
throughout, so we set density p — 1, and V ■ v = 0, where v is velocity. We include the 
thermal pressure P, and the magnetic pressure in the total pressure p = P + b 2 /2 and 
assume isotropic kinetic and magnetic viscosities, v v . The MHD equations become, 



d t v = -v • Vv - Vp + b • Vb, +u v V 2 v, 



(1) 
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d t b = -v ■ Vb + b ■ Vv + u b V 2 b, (2) 

To relate p to v and b, we take the divergence of equation [1], which upon inversion, yields 

r d 3 x' (Vv : Vv - Vb : Vb) 

P = J 17 • (3) 

A random forcing field with energy 6f ■ /St is generated at each time step and added 
to the existing velocity field, where e/ is the average forcing power. The amplitudes of 
the forcing field Fourier modes are assigned according to a specified power spectrum, with 
the energy selected from a Boltzmann distribution. The mode phases are random within 
the constraint of divergencelessness. We input kinetic helicity v ■ V x v at the forcing 
wavenumber of s = k/27c = 4.5 by making a randomly determined subset of the Fourier 
modes maximally helical, leaving the rest unchanged. The fraction of maximally helical 
modes is fh, which we denote "fractional helicity." In contrast, simulations of Maron & 
Cowley (2001) invoked zero mean magnetic field and zero mean kinetic helicity. Only 
fractional random fluctuations of the kinetic helicity of order 10% were present. The 
magnetic helicity was also initially zero and subsequently fluctuated about zero at an 
amplitude of 10% of the maximum. The equations of MHD are solved spectrally. The 
turbulence is incompressible and the boundaries are periodic. (Including compressibility, 
while ultimately important for detailed applications, is not expected to have a dramatic 
effect on the qualitative conclusions herein regarding the role of helicity.) Wave numbers 
and physical scales are related by Xk = 2tt. Viscosity and resistivity are of the k 2 type 
{i/ v V 2 v and z/f,V 2 6). The code is exhaustively discussed in Maron & Goldreich (2001). 
We note that at each time step, the time derivative of the magnetic helicity is equal to 
the volume integral of the current helicity times the resisitivity. The helicity conservation 
equation is satisfied and this is important in what follows. 

The other key parameters are as follows: the magnetic Prandtl number is 
Pr = v v /ub ~ A^/A^ , the ratio of the viscosity to magnetic diffusivity, where X Uv and X Ub , 
are the viscous and resistive scales respectively. We denote v\ and b\ as the speed and 
magnetic field at scale A, and Vf and Xf as the forcing scale RMS velocity and forcing scale 
respectively. When Pr > 1, the scales have the ordering A/ > X Vv > X Ub . 



3. Results 



We show results here for a selection of 64 3 simulations which is sufficient to identify 
the basic effects of fractional helicity on the location of energy peaks. The fh ranges from 
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to 1 by increments of 0.1 for simulations AO through A10. For all A0-A10, we used a 64 3 
grid, Sf = 4.5, u v = 3 x 10~ 3 , v h — 1 x 10 -3 , and Pr = 3. 

The usual kinetic and magnetic energy spectra are defined as the quantities inside 
the energy integrals E v = J E v (s)ds and E b = J E b (s)ds respectively. The spectra for a 
range of values of fh are shown in Fig 1. The time evolution of the fh = l case is shown in 
Fig 2 and the time growth of magnetic helicity is shown in Fig 3. Notice in these figures 
that for fh > fh,crit ~ 0.5 the peak at the forcing scale grows as does the large scale field. 
For fh < fh,crit, the large scale field decays, no peak appears at the forcing scale, and the 
magnetic helicity in the box grows very weakly, if at all. Though we have presented only 
Pr = 3 cases in the figures, we have also performed simulations with Pr = 9 and found 
that fh,crit ~ 0.7. Thus we find that fh, C rit increases with Pr. 

We checked for hysteresis by using the saturated state of the fh — 1 simulation of Fig. 
3 as the initial condition for another simulation with fh = 0.4. We found that the magnetic 
helicity subsequently decayed to the same value as in the simulation which started with a 
weak mean field with initial fh = 0.4. There was no evidence for hysteresis. 



E„(k), f h =0 

E b( k ). f h=- 4 
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Figure 1: Saturated kinetic and magnetic energy spectra for successive values of fh- The 
simulations are AO, A\, A5, and A10. The kinetic spectra are identical for f h < 0.4 
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Figure 3: Evolution of the magnetic helicity as a function of fh for simulations from series 
A. At t = 7, simulation A10 (fh — 1) has developed a large scale magnetic field. 
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4. Discussion 

Forced helical turbulence was studied systematically in Pouquet et al. (1976) using the 
Eddy Damped Quasi- Normal Markovian (EDQNM) closure scheme (Orsag 1970). Forcing 
with kinetic helicity leads to segregation of magnetic helicity because the magnetic helicity 
growth equation has a source term that depends on the kinetic helicity The kinetic helicity 
inputs one sign of magnetic helicity at small scales but the opposite sign is generated on 
large scales. The growth of the large scale field, represents a non-local inverse cascade of the 
magnetic helicity from the forcing scale and can be interpreted as an a 2 MFD (Pouquet et 
al. 1976; Brandenburg 2001; Blackman k Field 2001b; Field k Blackman 2002; also Verma 
2001). Averaged over a periodic box, the time evolution of the total magnetic helicity 
satisfies <9j(A ■ B) = —2ub(J ■ B), where A is the vector potential and J is the current 
density. If we divide (A ■ B) into small scale and large scale contributions, we can see that 
after the large scale helical field energy grows to k s B 2 /k[, where s,l refer to the dominant 
small and large scale, the large scale magnetic helicity dominates (Fig. 3). The growth 
saturates as long as there remains a net current helicity, but the growth rate is resistively 
limited, implying a "slow" (decreases with increasing k Vb ) MFD. The large scale field growth 
of Fig. 3 for all f\ > fh,crit is sl° w i n this sense. (Note however that there is a short initial 
"fast" phase. Field k Blackman (2002) and Blackman (2002), discuss this in the context 
of a dynamical quenching model which fits Brandenburg (2001) and the results herein, and 
the asymptotic quenching of Gruzinov k Diamond 1994; Bhattacharjee k Yuan 1995). 

The sign of the magnetic helicity of the growing large scale field is opposite to that of 
the kinetic helicity. This is consistent with MFD theory if the kinetic helicity dominates the 
a effect of the MFD: A positive kinetic helicity means that a would be negative. But the 
growth of the magnetic helicity associated with the large scale field is oc aB 2 (Brandenburg 
2001, Blackman k Field 2001; Field k Blackman 2002) so that a positive input kinetic 
helicity, which gives a negative a, produces a negative large scale magnetic helicity. Fig. 3 
shows that (A • B) is dominated by the large scale contribution. 

That the large scale field growth proceeds asymptotically "slow", would seem to 
threaten the relevance for e.g. the Galaxy as it is commonly argued that the MFD of the 
Galaxy has to be "fast." (c.f. Ruzmaikin et al 1988; Zweibel k Heiles 1997). Unlike a 
periodic box, real astrophysical rotators have boundaries, shear, and have helicity driven 
by the combination of the underlying rotation and stratification with a spatial variation 
in transport coefficients. So these differences will ultimately need to be studied before 
results from periodic box solutions like ours can be directly applied to astrophysical large 
scale field growth (Blackman k Field 2001). Also, the large scale field for > fh,crit i n 
the simulation becomes super-equipartition as seen in Figs. 1 & 2, because it is nearly 



-7- 



force-free. That being said, the enterprise of periodic box simulations is very worthwhile; 
even comparing the difference between periodic boxes with more realistic simulations will 
be important to our further understanding MHD turbulence and allows the demonstration 
of some principles that emerge in the simplest possible forced non-linear dynamo system. 

In spite of these caveats, the results of the present simulations are provocative: the 
same fh,crit determines both whether the large scale field grows and whether a peak grows 
at the forcing scale. The large scale field growth and the presence of the small scale peak at 
the forcing scale appear to be intimately related. Magnetic helicity undoubtedly also plays a 
role drain of the small scale peak from the resistive scales for fh > fh,crit- Consider the role 
of magnetic helicity conservation in mode interactions via a slightly more general version of 
the argument of (Frisch et al. 1975): Suppose an initial state of maximal magnetic helicity 
is confined to wave numbers k m and k n , with k m < k n and suppose the magnetic field 
dominates the energy at these wavenumbers. We then have E b (k m ) + E b (k n ) = E T (k p ) and 
Hb(k m ) + Hb(k n ) = Eb(k m ) I k m + Eb(k n )/k n < ET(k p )/k p , where Et means the total energy. 
The last inequality can only be satisfied if k p < k n . This argument applies only when the 
scales k m and k n are magnetically dominated, though k p need not be (a point not addressed 
in Frisch et al. 1975). An initial state for which the field is dominant, and therefore satisfies 
the validity criterion, is the small scale saturated state shown in Fig. 1 for fh — 0. We can 
reason that when sufficient magnetic helicity is imposed, it, and its associated energy would 
drain from the small scales, at least until the field reaches equipartition with the velocity. 
This is qualitatively consistent with the observed deficit in the magnetic energy from the 
large k region in figures 1 and 2 for f h > fh, cr it, as compared to f h < f h , C rit- 

The growth of the actual peak in magnetic energy at the forcing scale is also aided 
by the fact that a forced kinetic helicity reduces the non-linear transfer term in the 
Navier-Stokes equation. The non-linear term — v • Vv = v x u — Vv 2 . When helicity 
is present, the v x u term is reduced. For sub-sonic turbulence, the v 2 contribution to 
the evolution equation should be inconsequential. Thus, since the main cascade driver is 
reduced for helical turbulence, the kinetic energy requires more time to cascade, providing 
a bit more time for this energy to be transfered directly into stretching the magnetic field 
near the input scale. Though a cascade of magnetic energy steadily drains the field from 
the input scale, the hold up of the kinetic energy cascade means that there is more time to 
resupply the field to a larger amplitude before draining, compared to the fh < fh,crit case. 

The total kinetic energy density for k > kf is fixed and is always larger than the total 
magnetic energy density for k > kf. If fh < fh,crit, then there is a significant non-helical 
part of the field which feels no tendency to inverse cascade. This fraction piles up quickly 
on the small scales (Chou 2001; Maron & Cowley 2001). If this fraction dominates, then 
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the spectrum will be dominated by the non-helical turbulence dynamics. If fh > fh,crit, 
the magnetic energy associated with the helicity, which inverse cascades, dominates. A 
lower limit on fh,crit can De found from modeling the large scale field growth as an a 2 
dynamo, which has a growth rate of a — /3s\, where f3 is turbulent diffusion and s\ — 1 is 
the growing large scale wave number. But initially a and (3 have their kinematic values and 
cx/si/3 ~ 2f h Sf/(3si), and so initial growth requires at least fh > si/sf ~ 0.33 = f crit . This 
is roughly consistent with our results within small factors of order 1. 

At risk of applying our "slow" dynamo spectral shape results too cavalierly, we consider 
the implications for the Galaxy. The Galactic field has both a large scale (> 2kpc) and a 
small scale component (< lOOpc) (e.g. Zweibel & Heiles 1997). The small scale field (which 
has magnitude ~ few times larger than the large scale field) appears to have a peak at the 
forcing scale, as does the kinetic energy (Armstrong et al. 1995), and the two are in near 
equipartition with v ~ b ~ lOkm/s. Ignoring the "slow" vs. "fast" issue for the moment, 
our results would imply that the Galactic magnetic spectrum (with its small scale peak 
at the forcing scale) is only consistent with turbulence forced with f\ > fh,crit- Note that 
we have been considering only externally forced turbulence, as opposed to self-generated 
turbulence from shear. This particular assumption, at least, is consistent with the Galaxy, 
where supernovae are the primary driver (Sellwood & Balbus 1999). (Given that our box 
dynamo is "slow" we should note that well motivated analytic interpretations for generation 
of "fast" large scale magnetic energy in sheared rotators which appeal less explicitly to 
kinetic helicity, or not at all have been studied (see Balbus & Hawley 1998; Vishniac & 
Cho 2000), as do specific proposals for how boundaries might enable fast helical a — VL 
dynamos (Blackman 2002). We do not discuss these further but note that in the Galaxy, a 
net magnetic flux in the Galactic disk in addition to magnetic energy seems to be needed.) 

Our results might also suggest that kinetic helicity would play a role in the protogalactic 
small scale dynamo model (Kulsrud et al. 1997) of the large scale Galactic field. In this 
model, the large scale field of the Galaxy results from gravitational collapse and flux freezing 
of the small scale protogalactic field. The model requires that small scale dynamos generate 
significant power at the forcing scale, and our results would suggest this is only possible in 
Pr > 1 plasmas when f h > f h ,crit- 

5. Conclusions 

The magnetic spectrum of MHD turbulence forced in a periodic box with fractional 
kinetic helicity above a critical value fh, C nt) saturates with two peaks: a large scale peak 
and a peak at the forcing scale, when Pr > 1. If fh < fh,crit there is only one peak 
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in the spectrum, and it is at the resistive scale. Though the turbulence is not strictly 
Alfvenic for any f^, it is much more nearly Alfvenic for larger f^. The range of scales, the 
boundary conditions, and the actual nature of helical forcing pose obstacles when comparing 
non-sheared periodic box simulations with real astrophysical rotators. Nevertheless, the 
important qualitative implication of our results is that helical forcing not only influences 
the large scale magnetic field spectra of forced turbulent systems but may also help account 
for observed peaks at the forcing scale. 
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